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Abstract: The activation probe is a robust tool to measure flux of fusion products from a magnet¬ 
ically confined plasma. A carefully chosen solid sample is exposed to the flux, and the impinging 
ions transmute the material making it radioactive. Ultra-low level gamma-ray spectroscopy is used 
post mortem to measure the activity and, thus, the number of fusion products. 

This contribution presents the numerical analysis of the first measurement in the ASDEX Upgrade 
tokamak, which was also the first experiment to measure a single discharge. The ASCOT suite of 
codes was used to perform adjoint/reverse Monte Carlo calculations of the fusion products. The 
analysis facilitates, for the first time, a comparison of numerical and experimental values for abso¬ 
lutely calibrated flux. The results agree to within a factor of about two, which can be considered 
a quite good result considering the fact that all features of the plasma cannot be accounted in the 
simulations. 

Also an alternative to the present probe orientation was studied. The results suggest that a better 
optimized orientation could measure the flux from a significantly larger part of the plasma. 
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1. Introduction 

The fast particles, such as fusion alphas, ICRH, as well as NBI generated ions, must be confined 
in magnetically confined fusion plasmas until they have released their energy. Therefore it is es¬ 
sential to measure and understand the fast particle behaviour in burning plasmas, such as in ITER. 
Detecting the escaping particles is the first step in learning to confine them. 

The fusion product activation probe is a member of the family of probes and detectors that 
measure the flux of charged fast particles near the first wall of a magnetic confinement fusion 
device [HBB. The fast ions are often measured with scintillation materials, Faraday cups and 
semiconductor detectors which can measure on-line the flux with high temporal resolution. In 
contrast, the activation probe is based on post-mortem measurement of the irradiation induced 
activity of solid material samples. The flux of fast particles transmute part of the target material 
making it radioactive. The resulting gamma activity can be measured with ultra-low level gamma- 
ray spectroscopy. The measurement is well calibrated, which makes the probe particularly well 
suited for validation of fast particle codes and for cross calibration of other detectors in present 
day tokamaks. The activation probe is very robust, which make it uniquely suitable for the harsh 
environment of a fusion reactor. The activation method requires neither optics nor electronics inside 
the plasma chamber, where the neutron and gamma radiation would inevitably deteriorate them. 
The probe is also inherently isotope selective and practically insensitive to the neutron flux. The 
measurement can be inherently energy selective, when using an activation process with a threshold 
energy. 

The fusion product flux at the midplane manipulator in ASDEX Upgrade (AUG) was measured 
for the first time [0] using an activation probe [0]. This contribution presents numerical analysis of 
these experiments using the ASCOT code [0]. In this work, the fusion reactivity is calculated using 
ASCOT and its AFSI Fusion Source Integrator component, and for verification purposes, also with 
TRANSP [0]. The detailed probe geometry was included in the AUG 3D first wall model. We 
performed an adjoint Monte Carlo simulation, that was equivalent to following the collisionless 
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Figure 1 . The poloidal cross section of the plasma. The flux surfaces are indicated in green and the activation 
probe in red. 

orbits of the fusion products from their birth location in the plasma to the probe. In the adjoint 
method, markers are launched backward in time from the probe towards the plasma. Only fusion 
protons were experimentally measured, but tritons and helium-3 were modelled for comparison. 

The main result of this paper is the comparison of measurements and simulations to demon¬ 
strate quantitative analysis capability. Indeed, the current analysis is in reasonable quantitative 
agreement with the measurements. The presented numerical analysis provides, for the first time, 
absolutely calibrated flux of fusion products to the probe. This is an improvement compared to the 
previous analysis [0 performed on JET measurements, that only provided relative profiles of the 
flux as a function of major radius and thus required external calibration. Furthermore, the analysis 
was performed for an alternative probe geometry, as suggested by a previous study 0. 

2. The probe, the measured plasma and the numerical methods 

An activation probe exposes carefully chosen material samples to the escaping fusion product flux 
from the plasma. In the AUG xperiments, the probe was connected to the midplane manipula¬ 
tor slightly above the midplane (figure [TJ). The geometry is similar to that foreseen for an ITER 
activation probe [01- 

In the experiment, the probe was exposed to the dedicated ASDEX Upgrade discharge #29226. 
The plasma flat top was stationary and approximately 6 s in length (figure @). It was a deuterium 
plasma with 7.25 MW of injected deuterium neutral beam power. The core temperature was ap¬ 
proximately 3keV and density 7 • 10 19 /m 3 . The plasma profiles were acquired using integrated 
data analysis UTIll for the time window 2.1-5.0 s, while the magnetic background was calculated 
for the time point 2.14 s. In the modelling, nitrogen was used as the effective impurity, to match 
the measured effective charge number. 
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Figure 2. Time evolution of the plasma parameters. The top panel shows the long, uniform plasma flat 
top. The second panel shows the neutron rate, both measured and one obtained from a transport calculation 
with TRANSP [Of The third panel shows how the neutral beams dominate the heating of the plasma. The 
bottom panel displays the line averaged electron density from interferometry channel passing through the 
core plasma and ion temperature from charge exchange recombination spectrometry. The shaded area is the 
time window for averaging of the background profiles from integrated data analysis 111 111 . The gray vertical 
lines mark the three TRANSP time points, that will be described in section 0. 


2.1 The activation probe design and operation principle 

The probe consists of a graphite sample holder (figure 0a) and a protective cover with a slit allowing 
fusion products to hit the samples (figure 0b). The slit was situated at the top of the graphite cover. 
Most of the hardware was reused collector probe components fUTH . Only the sample holder was 
manufactured for these experiments. The holder has four slots, spaced 90° apart, with room for 
six samples each. The holder can be rotated between discharges to allow measurement of multiple 
discharges. The sample materials were in these experiments boron carbide B 4 C, calcium fluoride 
CaF 2 and yttrium orthovanadate YVO 4 . The two B 4 C samples (second and fourth, when counting 
from the tip of the probe) are of interest in this study. The flux of protons to these samples from the 
D(D,T)p reaction was successfully measured, by measuring the amount of 7 Be isotope born in the 
l 0 B(p,o:) 7 Be nuclear reaction. It is this measurement that will be compared to the simulations. 

In the numerical analysis, the model of the probe was included in the detailed 3D wall (fig¬ 
ure 0c) of the tokamak. The general shape and dimensions of the graphite cap and sample holder 
with samples are included in the model as shown in figure 0 (d-e). 

2.2 The adjoint integration scheme 

The quantity of interest is the flux of particles from the large plasma (source) to the tiny samples 
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Figure 3. Illustration of the probe, (a) The sample holder with samples numbered, (b) The graphite cap 
covering the sample holder, (c) The 3D-wall model, including a model for the probe. The coloring indicates 
fusion product flux distribution onto the first wall as calculated by forward Monte Carlo ®|. (d) Close up 
of the probe geometry used in this study. The protective cover is shown in semi-transparent yellow. The 
uniquely coloured samples are partly visible through the slit (indicated by the arrow), (e) Cross-sectional 
view of the probe. The green mesh is used at the ends of the probe and the blue mesh is used at the ends of 
the cavity between the sample holder and the graphite cap. The samples and the opening in the cover are at 
the far right side. 


within the probe (target). The normal, forward in time, calculation is very inefficient, because only 
very few markers launched from the plasma will find their way to the target that is only visible 
through the narrow slit in the graphite shell. This can be remedied by swapping the roles of the 
source and the target in the adjoint Monte Carlo integration. In this approach, the markers start 
backward in time from the target and are likely to pass through the plasma before hitting a wall. 
In practise, the calculation is done by calculating two independent quantities and subsequently 
multiplying them to acquire the final flux. First the fusion product source rate is calculated and 
then markers are followed backward in time from the samples to the plasma. Both calculations are 
done assuming steady state. 

The calculation of fusion reactivity starts with the calculation of the NBI ionisation with the 
BBNBI code [Ql . The slowing down distribution of the guiding centers of the beam ions is sub¬ 
sequently calculated with ASCOT [0, taking into account collisions with the background plasma. 
The fusion reactivity is then be calculated using the plasma profiles and slowing down distributions 
using the AFSI Fusion Source Integrator, a part of the ASCOT suite of codes. The reactivity is 
stored in a ( R,z ) (major radius, elevation) grid. Using identical plasma profiles and slowing down 
distributions, a very good match (within a few %), was found between fusion reactivities calculated 
by AFSI and TRANSP [0] (figure @). 

In the second phase of the calculation, the markers are followed backwards in time (? prt ). The 
initial configuration space marker distribution Aw is uniform on the targets and the marker velocity 
distribution Q.\y is isotropic with a single velocity v. The full gyro motion of the particle orbits 
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Figure 4. A comparison of fusion reactivity calculated by AFSI and TRANSP/NUBEAM. The plasma 
profiles and fast ion distributions were identical in this comparison. Beam-beam (BB) reactivity corresponds 
to reactions that take place among beam ions, while beam-target (BT) reactivity reflects the fusions between 
thermal plasma and beam ions. Beam target rate is much larger than the BB rate in the studied plasma. The 
reaction shown here is D(D,n)He 3 . 

is followed backwards in time until they hit the probe structures or an other part of the tokamak 
wall. The adjoint density is the steady-state phase-space density of the reverse-in-time particles as 
integrated from their trajectories. The adjoint density is calculated in a similar grid as the reactivity. 
The flux of fusion products arriving at the target is then obtained by multiplying the two grids 
element by element and summing them up. This method approximates the following integral for 
the total flux: 



(2.1) 


fusion reactivity 


For more details, please see HBH . 

3. Simulation results and comparison to experiments 

The calculation of flux was done in three phases, calculation of the fusion reactivity, calculation of 
the adjoint density and multiplication of these two together. The fusion reactivity was calculated 
for thermonuclear, beam-target and beam-beam fusions. The sum of these three reactions is shown 
in figure ||(a) for the reaction D(D,p)T. The other reactivity, for the reaction D(D, ’Hc)n, is nearly 


identical. 


The adjoint calculation was performed by launching three million markers from the samples. 
Approximately 15% of the markers exited the graphite cap through the slit. The bulk of these mark¬ 
ers hit the wall in less than a microsecond, and produce the high adjoint density in the immediate 
vicinity of the probe, see figure |5](b). Even the longest living markers hit a wall within 100 mi¬ 
croseconds, during which they produce the volume of low adjoint density that extends to the centre 
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Figure 5. (a) The fusion reactivity, (b) The adjoint marker density of protons, summed up for all samples, 
(c) The origin of the contributing particles, (d) The adjoint marker density of protons if the probe were 
rotated by 90°. (e) The birth locations of the contributing particles. (90° rot.) 


of the plasma. The adjoint density is closely related to the “instrument function” of the probe: 
It directly indicates which parts of the plasma the probe measures and with what kind of relative 
sensitivity. In absolute sense, it indicates the “probing power” in phase space, which explains the 
somewhat counter-intuitively units, m 3 -srad. 

The total flux to targets was obtained by multiplying the fusion reactivity by the adjoint density 
(figure 0c), (and scaling with the relevant factors). The fusion reactivity drops quickly outside the 
core, while the adjoint density is large only near the probe. Consequently their product is non¬ 
vanishing mainly at the plasma periphery, at minor radius 0.7<p po i<l, at the low field side. As 
described above, the resulting flux is sum over the point-wise products and presented in figure 0(a). 
The figure also shows the measured flux, obtained by dividing the total number of particles by the 
flat top length of 6.0 s. Also a transport simulation of the discharge was performed with TRANSP 
to verify the forward part of the calculation and to assess the simulation errors. The calculation 
was repeated for three time steps from the TRANSP simulation (marked in figure @). The four 
calculations are within approximately a factor of two of the experimental measurement as figure 0 
reveals. 

Previous, forward Monte Carlo study [0] suggests that if the probe were rotated by 90°, the 
measured fusion product flux would be much higher, because many more products hit the sides of 
the probe than the top (figure [7]a). This was studied with the adjoint method, and indeed, rotating the 
probe by 90° strongly changed the instrument function (figure 0d). Rotation of the probe changed 
the sampled velocity space. In the rotated configuration, the detected products had more parallel 
velocity when hitting the probe. Such particles were able to evade the walls near the probe and 
enter from the core of the plasma. Therefore, the probe sampled the fusion reactivity not only 
at the edge, but also in the core, especially from the high field side 0.3<p po i<0.8. The plasma 
edge contribution from near the probe was reduced, but not completely eliminated (figure 0e). The 
combined effect was increase of the signal, but only in certain samples that received the flux from 
the core (figure [7|b). 
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Figure 6. (a) Simulated and measured flux to the samples. In each group of bars, the left one corresponds 
to 3 He, the middle one to protons and the right one to tritons. The red bar is the measured proton flux with 
uncertainty, (b) Comparison of measurements and adjoint simulations using four different fusion reactivities. 
The green area indicates the measurement uncertainty. 
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Figure 7. (a) Forward Monte Carlo simulation results |@| indicate that rotating the probe by 90° would have 
dramatically increased the flux to the samples. The left panel shows model of the probe graphite cover. Each 
circle depicts a fusion proton hit. The blue ones are approximately the ones that would enter through the 
slit (at the top of the probe) and green ones have hit the end of the probe. The right panel shows the side 
of the graphite cover unwrapped onto a plane. The maximum density of protons is situated approximately 
90° from the slit, (b) Adjoint Monte Carlo simulated flux to the samples in case of rotating the probe 90° 
clockwise as viewed from the plasma, resulting in a horizontal slit. In each group of bars, the left one is 3 He, 
the middle one protons and the right one tritons. 


4. Summary and Discussion 

We have demonstrated that the activation probe measurements at AUG can be quantitatively re¬ 
produced using a novel combination of regular Monte Carlo and adjoint Monte Carlo methods, 
realized by the ASCOT code in this contribution. A successful comparison between the absolutely 
calibrated measured flux and the numerically computed value was achieved for the first time. The 
modelling could be improved, by including some of the following features neglected here: the 
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simulations were performed for single time-slice while the plasma evolves in time. The fusion 
products were assumed to be collisionless. Their birth distribution was assumed to be isotropic and 
monoenergetic, which is certainly not true for beam-target fusion products. Regardless of these 
simplifications, the calculated fusion product flux to the probe was within approximately a factor 
of two of the experimentally measured uncertainty. 

The analysis of rotated probe head suggests that there is large room for optimising the mea¬ 
surement geometry and this contribution describes the tool for the studies. One could conceivably 
design a probe that is sensitive only to a specific part of plasma. It would be interesting to study 
the effect of pushing the probe nearer the plasma, rotating the probe, changing the slit shape etc. 
These studies could also be done for ITER, for which the robustness of the activation probe offers 
unique possibilities. 

Acknowledgments 

This work was partially funded by the Academy of Finland project No. 259675. This work has 
received funding from Tekes - the Finnish Funding Agency for Innovation under the FinnFusion 
Consortium. 

References 

[1] N. Gorelenkov, S. Pinches, and K. Toi, Energetic particle physics infusion research in preparation for 
burning plasma experiments, Nuclear Fusion 54 (2014), no. 12 125001. 

[2] M. Cecconello, O. M. Jones, W. U. Boeglin, R. V. Perez, D. S. Darrow, I. Klimek, S. E. Sharapov, 

M. Fitzgerald, K. G. McClements, D. L. Keeling, S. Y. Allan, C. A. Michael, R. J. Akers, N. J. 
Conway, R. Scannell, M. Tumyanskiy, G. Ericsson, and the MAST Team, Energetic ion behaviour in 
mast. Plasma Physics and Controlled Fusion 57 (2015), no. 1 014006. 

[3] G. J. Sadler, S. W. Conroy, O. N. Jarvis, P. van Belle, J. M. Adams, and M. A. Hone, Investigations of 
fast-particle behavior in joint european torus plasmas using nuclear techniques. Fusion Science and 
Technology 18 (1990), no. 4 556-572. 

[4] G. Bonheure, M. Hult, A. Fenyvesi, S. Akaslompolo, D. Carralero, D.Degering, A. de-vismess Ott, 

M. Garcfa-Munoz, B. Gmeiner, A. Herrmann, M. Laubenstein, G. Lutter, H. Mueller, V. Rohde, 

W. Suttrop, G. Tardini, and the ASDEX Upgrade team. First escaping fast ion measurements in 
ITER-like geometry using an activation probe, in Proc. 40th EPS Conference on Plasma Physics, 
Europhysics Conference Abstracts, (Espoo, Finland), p. 06.510, 2013. 

[5] G. Bonheure, M. Hult, R. G. de Orduna, D. Arnold, H. Dombrowski, M. Laubenstein, E. Wieslander, 
T. Vidmar, P. Vermaercke, C. P. V. Thun, M. Reich, S. Jachmich, A. Murari, S. Popovichev, J. Mlynar, 
A. Salmi, O. Asunta, M. Garcia-Munoz, S. Pinches, R. Koslowski, S. K. Nielsen, and JET-EFDA 
Contributors , Experimental investigation of the confinement of D(^He,p)a and D(D,p)T fusion 
reaction products in JET, Nuclear Fusion 52 (2012), no. 8 083004. 

[6] E. Hirvijoki, O. Asunta, T. Koskela, T. Kurki-Suonio, J. Miettunen, S. Sipila, A. Snicker, and 
S. Akaslompolo, ASCOT: Solving the kinetic equation of minority particle species in tokamak 
plasmas. Computer Physics Communications 185 (2014), no. 4 1310-1321, [jarXiv: 1308. 1904Q . 




[7] R. Budny, M. Bell, A. Janos, D. Jassby, L. Johnson, D. Mansfield, D. McCune, M. Redi, J. Schivell, 
G. Taylor, T. Terpstra, M. Zarnstorff, and S. Zweben, Simulations of alpha parameters in a TFTR DT 
supershot with high fusion power. Nuclear Fusion 35 (1995), no. 12 1497. 

[8] Akaslompolo, G. Bonheure, G. Tardini, M. Garcfa-Munoz, T. Kurki-Suonio, and the ASDEX 
Upgrade team, ASCOT simulation of fusion product activation probe experiment in ASDEX Upgrade 
tokamak , in Proc. 40th EPS Conference on Plasma Physics, Europhysics Conference Abstracts, 
(Espoo, Finland), p. P5.101, 2013. 

[9] G. Bonheure, M. Hult, R. G. de Orduna, P. Vermaercke, A. Murari, S. Popovichev, and J. Mlynar, 
Fusion alpha loss diagnostic for ITER using activation technique. Fusion Engineering and Design 86 
(2011), no. 6-8 1298 - 1301. Proceedings of the 26th Symposium of Fusion Technology (SOFT-26). 

[10] W. Schustereder, A. Herrmann, and V. Rohde, Impurity survey after boronization in asdex upgrade 
measured by a collector probe, Physica Scripta 2007 (2007), no. T128 14. 

[11] S. K. Rathgeber, R. Fischer, S. Fietz, J. Hobirk, A. Kallenbach, H. Meister, T. Piitterich, F. Ryter, 

G. Tardini, E. Wolfrum, and the ASDEX Upgrade Team, Estimation of profiles of the effective ion 
charge at ASDEX Upgrade with Integrated Data Analysis, Plasma Physics and Controlled Fusion 52 
(2010), no. 9 095008. 

[12] G. Tardini, C. Hohbauer, R. Fischer, R. Neu, and the ASDEX Upgrade Team, Simulation of the 
neutron rate in ASDEX Upgrade H-mode discharges. Nuclear Fusion 53 (2013), no. 6 063027. 

[13] O. Asunta, J. Govenius, R. Budny, M. Gorelenkova, G. Tardini, T. Kurki-Suonio, A. Salmi, S. Sipila, 
the ASDEX Upgrade Team, and the JET EFDA Contributors, Modelling neutral beams infusion 
devices: Beamlet-based model for fast particle simulations. Computer Physics Communications 188 
(2015) 33-46. 

[14] S. Akaslompolo, Adjoint monte carlo calculation of charged plasma particle flux to wall, April, 2015. 
larXiv: 1504, 00214| . 


- 9 - 




